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^^ ' We study two-leg spin ladder systems analytically by using a. t-J model which includes the in- 

O^ ' terchain spin exchange and the interchain hopping term. The spin part is mapped to a quantum 

(-H , sine-Gordon model via a bosonization method that fully accounts for the phase fluctuation. The 

j^ ' spin gap in the Luther-Emery phase evolves to gapless Luttinger liquid phases in even-ladders at 

certain doping concentration. This transition occurs as doping rate increases and/or the ratio of the 

interchain spin exchange and hopping to the intrachain couplings decreases. We also estimate the 

transition temperature at which the conventional electron phase can be deconfined to spinon and 

holon. 
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PACS numbers; 75.10.Jm, 71.27.-fa, 71.10.Hf, 74.20.Mn 



I. INTRODUCTION 



P^ . Low-dimensional quantum magnetic systems, such as spin chains, spin ladders, and planar lattices provide very 

O useful information about the physics of high Tc superconductors. During the past few years, spin-ladder systems 

, O, , have been studied extensively by various groups ||l| H] [3, especially by numerical methods or by mean-field (MF) 

calculations M M. A "spin ladder" consists of parallel "chains" of magnetically bound Cu atoms connected by 

psj \ interchain couplings ("rungs"), as found in 5'r„_iCu„+i02n, -^a4-(-4nC'M8-i-2nOi4+8n (Fig. 1). Spin ladder compounds 

^ • may become superconductors if doped with charge carriers [pj M. Recently, superconductivity has been found in 

CN I Sri4-xCaxCu24:Oii below 12 K and with external pressure of 3 GPa. 

^^ . Beginning with the half-filled ground state of a spin liquid, short-range singlets are formed predominantly across 

^-J ' the rungs of the ladders from the antiferromagnetic (AF) states. The unusual spin-liquid nature of an undoped parent 
j-^ , system and the evolution of the finite gap upon doping have also been studied by numerical methods |lQ] . Excitations 
^>^ • of the singlet ground state yield spin triplets which propagate along the ladder direction. The study of even-leg ladder 
^\ , systems has shown that spin gaps decrease as the number of legs increase. ( The gap remains finite down to J± = J\\ , 
"^ ' which is a signature of superconductivity pT| ] ) . Within the mean-field theory, in numerical works and renormalization 
Ci ■ group (RG) calculations, a superconducting phase with a d-wave nature emerges from the short-range RVB liquid 
2 . upon doping. Then, new quasiparticle (QP) states appear which carry both charge (positive) and spin. 

In the weak-coupling approach, one starts with two non-interacting spin chains and turns on perturbative interchain 
interactions. For finite interchain interactions, this method does not accurately describe the correlations near the phase 
transitions. Zhang et al. 112] have studied two-leg t ~ J ladders without the interchain exchange by mapping the spin 
Q part to a sine-Gordon model. Our study includes the spin exchange interaction across the rungs of the spin ladders 
and is based on a (nearly) exactly solvable model. Based on the soliton description, we discuss the phases and possible 
transitions between them. Exactly solving for spin ladders gives better understanding of the correlations between 
phase excitations and topological defects, while comparing with some weak-coupling studies Q] and MF results which 
^ ' partly addressed this issue 13 . 

We use a t-J model for two-leg ladders in the strongly(isotropically) coupled regime in which the interchain interac- 
tions {t±, J±) are larger(equal) than the intrachain interactions (iy, Jy). Introducing weak coupling between two-leg 
ladders is staightforward for even-leg ladder systems. Recent studies pSi have shown that ladder systems with an 
odd number of legs are conducting Luttinger liquids in the odd-parity channel and spin liquids in the even parity 
channel jlj] . Odd- legged ladders can be mapped to a single-chain system while even-legged ladders are equivalent to 
a double-chain ladder. 



II. THE MODEL 



An effective model for the t-J Hamiltonian is obtained as follows 



^ = -^11 E E i^LCl.a + h.C.) 
1=1,2 <ij>,cr 

~t^ T. T. {c+^ci,.. + h.c.) 

Ll' = l,2 <ij>,rj 
1=1,2 <ij>,(7 
l,V = l,2 <ij>,<J 

+ MEE^^tcz.. (1) 

1 = 1,2 <T 

where < ij > is for nearest-neighbor( j = i + 1) interactions, / is the leg index and a is the spin (||) index. In strongly 
coupled t-J models with a large Coulomb repulsion (large U/t and large Jl/Jii), we impose the single-occupancy 
constraint "^^ C^Ci„ =0 or 1 at each site. The electron chemical potential /i is added for the constraint in the 
form of a Lagrange multiplier. The charge(hole) concentration is related to /x(l — (5) = 1/(1 + e^^^). In order to 
study the spin part and the charge part separately, we write the electron operator as the product of a charge carrying 
holon hi and a pseudospin operator Si: Cii = PhiSfP^ = hiSf , Ci^ = PhiS~ P^ = hiS~ . The projection operator 
P = {l—Tii^-cr) reduces the four on-site states |0, t>, |0,i>, |1, t>and |1, J,> to the physical Hilbert spaces \hole >, | t> 
and I i>, thereby imposing single-occupancy at each site. Since the charge sum rule (^^ C^Cia = 1 — <5) holds when 
the projection operators are neglected, we equate the electron operator with h^ S. The pseudospin operators acting 
at site i turn into spinless fermion operators a via a Jordan- Wigner transformation. 

III. CONTINUUM LIMIT AND BOSONIZATION 

In the spin exchange part, the pseudospin operators anticommute on the same site and commute on the different 
sites, so, we have H{J) parts become Jy ^i^i{S^S^_^i + h.c.) and J± J^i li^u^u+i + h.c). From Jordan- Wigner 
transformation, the intrachain term S^ S^^-^^ — a^ exp[— i7ra^ai]ai+i = af {1 — 2af ai)ai+i for spin ^-system, while the 
interchain term S^iS~i_^-^ = a/',a_i^;+i — 2a^;a^;ai^;ai.;+i. We introduce spinor fields (pa for the slowly-varying field o^ = 
Oi/i" where a=even at even(2s) site i and odd at odd(2s+l) site i. The even-site and odd site are coupled by nearest- 
neighbor interaction, yielding ^^^a+S,+i = ^£o[0+(2s){<?!)o(2s-fl)-(?!>o(2s-l)}-f(/)+(2s-Hl){</)e(2s)-(/)e(2s-2)}]. 
In continuum, the fermion field is given 4>a{s) — \/2ao^ a{x) where x = 2sao with the lattice constant oq. Then, 
(f>o{2s + 1) — (j)o{2s — 1) w 2aoy/2aodx'i'o{x)- The commutation relations are {^^(x), ^/3(2/)} — SapS{x — y) and 
{4>a{n), 4>p{ni)} = SapSnni- With the relation for exponentiated operators e^e^ = e'-^+^l (i-'[A.B]/2 ^ ^b ^a ^-[a,b] ^ 
we map the spin exchange interaction terms 

J\\ / ^ SiijJja ^ •^ll I'-'i S.j + "Jj Oj -\- OizOjz) 
<ij>,a 

= J||(a+ai+i - 2afa^a.iai+i + h.c.) = H{J\\). (2) 

Now the intrachain spin exchange part of Hamiltonian becomes 

H{J\\) = 2aoJ|| I dx 

[ ■^+{x)ia^d.,^,-,{x) -f *+(x)*^(a;)*+(x)z(Tia,*^(a;)] (3) 

where ^^i is the Fermi field of two species(bonding, anti-bonding band) and a^ is the Pauli matrix. 

Then, we bosonize the Fermi field with phase fluctuation parameters which can vary depending on the interactions 
(rather than the parameter being a constant VItt) . Employing the operator expansion of the Fermi field |1^ , \1/^ (cc) = 
eyipliaipnix)] exp[i/3/_ dydoipniyo^ y)] where the Bose field and its canonical conjugate momentum 11 = OqcJ) satisfy 



[(/)(x), n(x')] = i6(x — x') and [n(a;), 9o(/>(a;')] = -dx5{x — x'). In order to satisfy the commutation relation of ^ri{x), 
the phase fluctuation parameters a and /3 are related as 2af3 = ±7r. The Fourier components of density operators 
Pvil) = Z]/c*^(fc + '7)*»?(fc) satisfy the boson commutation relation [pRiq), PR{-qf)] = ^q 6qq,, [pL{q), PL{-q')] = 
■^q Sqqi and [pi?(g), PL{-q')] = 0, where prt^{x) ~ PLr){x) = H^ and pr,^[x) + pLr,ix) = ^d^cp^j is the fermion density 
deviation from its average. Introducing the dual field /_ dyllri{y) = -9rj{x) and linearizing the four branches of the 
spectrum near the Fermi surface(±fci?), "^nix) = e^'^^^Rri{x) + e~^^''^Lr^{x) with the right-going (R) and left-going 
(L) components for each specie, 

^vi^) = r?y P^ exp[-i(6'^(a;) - (?!)^(a;))], 

L^{x) = ^-— F^ exphi(6'^(a;) + </>^(x))] (4) 

where the Klein factors -F^ = exp[— z0^] are ladder operators [|l6| [ |l7| , raising/lowering the 77- fermion number by one. 
They are taken to yield the same signs for the species(?7) thus ignored in the thermodynamic limit, since the transverse 
hopping and spin exchange lead to interchain "pair fluctuations y] with both specie participating. Also, the '^ fields 
are written as even ^e?? = ~^(^>) + ^»?) ^-nd odd ^o») = ~7s(^-R)) + ^»))- Resultantly, the intrachain exchange J|| maps 
to the form of (9^0^)^ -I- 2cos(|0^). 

From the quartic interchain exchange J \_ part of Hamiltonian, —2J±a^ia^iai.iai^i+i = 

There are actually ten terms giving non-zero contribution to the resultant sine-Gordon, three of which arc shown 
below. 

H{Ja_) = -Ji_ ^ Su^S2ia = -2Jj_aQ / dx 

\ L+dxR2 - RtLiR+dxL2 + R2+L2L+d:,Li + ..]. (6) 

The terms in the integrand above break the left-right (chiral) symmetry which arises from the Umklapp scatter- 
ing between different chains (and backscattering between opposite spins). The interchain hopping term H{tj_) = 
—t± '^i„{Ci^^C2iiT + H.C.) yields identical forms in which we put the intrachain charge hopping OP /i; j/i;,i+i = a(5 = 

(5|| where < a < 1 and interchain OP ^; j/i;+i.i = (5[1 + 7^(1 — a)] = S±. Actual values of a, d±, Sn and above rates 
depend on the interaction strengths tj_, <||. 

Thus, the spin part of the Hamiltonian collected from the intrachain and the interchain contributions is 

Hs = Jdx ^[ (dl^i + dl4>2) + Wcos{^){i, + ^2)] (7) 

where (^ is a rescaled field from (f) = {1 + 2Wtt)(P with the lattice constant oq = 1. This is a quantum sine-Gordon 
equation which is the equation of motion of solitons (which are magnons with spin 1 (charge 0) made of two spinous 
bound by a confining potential). 

IV. PHASE FLUCTUATION AND THE INTERACTION STRENGTHS 

The constant prefactor is the ratio of interaction strengths interchain to intrachain 

^ \ t \ — J \ 

and is related to the phase fluctuation parameter /3^ = 167r/(l + 2W). In case of frustrated AF (which may occur in 
the ground state of odd-legged ladders), the sine-Gordon equation simplifles to H = -^ [(9^0)^ -I- W cos{(3(j))]. For 
any finite value of /3, the system has a zero-point fluctuation. In the dilute-soliton limit, the quantum fluctuation is 



dominated by the hard-core interaction between the spinless fermions 16 1. Depending on the parameter /3 or W, we 

J± 



have various soliton phases : (1) W — ^ = Wq : This is undoped state whose ground state is a AF spin liquid. 



as known from existing numerical studies ||l^. {2)W > Wq : Here we have Luther-Emery hquid with the gapped 
spinon(sohton) and gapless charge excitations. In the lower-temperature regime, excitonic soliton-antisoliton bound 
states exist. These correspond to integer spin magnon states which are bound states of two spin-i spinons. The CDW 
inherited from chains exists as well as the Luttinger liquid phase holons. When temperature increases, the excitonic 
bound states can be broken to spin-triplets (5 = 1 magnon)which make the system spin-gap metal. (3) ^ < M^ < Wq 
: In this case, repulsive interaction between spinons exists, and there are no stable solitons. This is equivalent to the 
strong fluctuation of a gauge field. SDW is in the background. (4) < M^ < ^ : The spinons are a Luttinger liquid 
as shown in algebraically decaying spin-spin correlation given below. The charge part is a Luttinger liquid as in (2). 
(5) W = or S±J± — t± : This limit corresponds to single soliton(spinon gas) mode. {Q)W = oo or 6\\J\\ = t|| : 
This corresponds to the classical limit with a finite soliton density. The electron states have dimensional crossovers 
from the incommensurate ID regime ior W < ^ to commensurate 2D-like states for Wq < W. At low doping and at 
low temperatures, the system can develop superconducting channel llG]. An excited spin phase generates magnons 
from triplets and solitons which propagate in the background of spin-density-wave(SDW). Weak but finite anisotropy 
or spontaneous breaking of chiral symmetry cause spin stiffness and low-energy excitations. The interplay of the 
excitation modes, such as order-preserving collective magnons and order-destroying topological solitons, may possibly 
lead to some long-range or intermediate-range order. Considering the asymptotic behaviours shown from scaling 
equations (Appendix), a schematic phase diagram is given in Fig. 2, which shows the various spinon and electron 
phases with the parameter W. 

V. CHARGE PART 

The separated charge part of the Hamiltonian is diagonalized via a Bogoliubov transformation of the holon operators. 
H = ^^.Ek{f^ fk -I- constants) where the holon excitation spectrum Ek = — 2i||Cos k^6j_tj_ is valid for tj_ > i||(l — 
cos 2STr) from the Luttinger sum rule. The charge phase with gapless excitation forms a Tomonaga-Luttinger ground 
state. According to numerical works fl^ [O, holes have lower energies when they form pairs across the rungs of 
the ladders. Thus, low-energy excitations include hole-hole pairs and hole-spin bound states (with positive charge) 
propagating along the ladder with a kinetic energy of about 4tn/ Jl. These excitations as well as the fiuctuations of hole 
pairs move against the lowest-energy background of charge-density wave(CDW). The holon speed is v^ — 2tnsinSTT. 

VI. CORRELATIONS 

The asymptotic electron Green's function can be calculated from perturbation theory excluding the vertex correc- 
tions ||l^ iQ and is given by 



where 7 — {y/l + W^ — l)/2 with a spinon velocity Vs — vp ^ I and ry = 1/2. The propagator Ge can be determined 
by convoluting the holon propagator and the spinon propagator and has no simple poles for low-energy states. Thus, 
the Fermi surfaces are considered separately for the holon and for the spinon. The size of the holon Fermi surface 
is found from J2k ui^hi^^^) — ^ (^ being the doping concentration). The ground state of the half-filled insulator 
is a disordered (or dimerized) spin liquid made of spin singlets. These singlets which form across the rungs become 
spin triplets as they are excited. Also, the triplets can propagate along the ladder as well as along the chain. The 
propagation across the rungs and ladders are "incoherent" in the sense that there is no stable eigenstate for the 
electron. The quantum numbers, such as spin, charge, and momentum, are not the same throughout the propagation 
from one leg to another [gl|. As separated phase fields, the charge and the spin degrees of freedom propagate with 
different speeds. At small lengths, the phase separation may manifest itself spatially in real space, for example as 
striped phases, as discussed by some authors ||2^. To see whether there is any strong-coupled fixed point, consider 
turning off the interchain hopping t±. The renormalization group (RG) dimension of the wavefunction comes from 
rescaling the variables {x,y,t) —>■ {ax,y,at). From the propagator, we have ip — > a~2-^. So, the Hamiltonian H± 
has the RG dimension of ^ -I- 77. The interchain hopping is irrelevant and the Luttinger liquid becomes a fixed point, 
if and only if r/ > 2 ■ This corresponds to the case of interchain interactions being much weaker than intrachain ones 
(t_L ^ i||). However, in strongly coupled ladders, Ja_/J\\ remains sufficiently large to develop gaps |g^ ||2J| and RG 
flow is away from the unstable Luttinger liquid fixed point B. Given that t± is greater than the spin gap, there is no 
strong-coupled fixed point, for even-legged ladders(?7 < ^). The equal-time spin-spin correlation is a power-law form 



< S{r)S{0) >w x-^^+^" cos 2kFX, (10) 

which shows the spinons being a non- Fermi hquid. From second-order perturbation theory |25[| , the spin gap size is 

J± + J\\cos k\\ + ^ J|?/,/l(3 — cos 2fc|[) with fcy = vr. 

VII. DECONFINEMENT 

In the 2D AF Heisenberg model or CP^ model, deconfinement of the electron (to a holon and spinon) occurs at 
T = |26|. As for spin-ladder systems, the "electron" phase is realized by having the spin and the charge phases 
confined at finite T above the transition temperature Ten- Since the confinement-deconfinement(CD) transition 
occurs in the absence of ordering, it is like a Kosterlitz-Thouless transition occuring in (2-f 1) dimensional system. 
From the partition function, Tcd is found 

< Tcd < TT^^ h (11) 



2k[6j_tj^ - J_l] '"^ k[Sj_tj_ - J_l] 

where k is Boltzmann constant. Hence, for T > Tcd, the electron propagator has the product form Gs{r,t) * 
Gh{—r,—t). Below Tcd, the spin phase and the charge phase are separated due to the screening of the confining 
potential between holons and spinons. Spin-charge separation is a well-accepted fact for ID systems which are 
Luttinger liquids. Spin ladders make transition, at some doping rates S, to a Luttinger liquid phase, where both spin 
and charge excitations become gapless. 

In summary, we map the spinons to solitons which are either magnons (excitons when bound) or excited phasons. 
From the soliton gap, we have the spin-gapped mode while the charge part is gapless Luttinger liquid. As the doping 
increases the ladders can develop superconducting phase or lose the spin gap beyond certain hole concentration. Also, 
depending on the ratio of interchain couplings to intraleg interactions, the spin ladders can become gapless Luttinger 
liquid. Our phase diagram is consistent with the features found in MF studies or weak-coupling approaches, such as 
presence of gapless mode and superconducting phase from even-ladders. The comprehensive diagram shows possible 
transitions as an explicit function of coupling ratio and doping. The electron propagator and the spin-spin correlation 
obey power-laws with the exponents depending on the interaction strengths and the doping rate. This shows the 
correlation of topological defects and phase excitations. Constructive interferences of phase excitations (solitons and 
magnons) may lead to the formation of intermediate-range or long range order, which is our forthcoming study. The 
electron is seen as a confined phase of the gauge flux associated with each point. Below the estimated transition 
temperature, the spin phase and the charge phase are separated. Henceforth, even-ladders show the crossover from 
gapped 2D-like states to gapless ID-nature, depending on the interleg couplings and/or the doping rate. 
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APPENDIX A: RENORMALIZATION GROUP EQUATIONS 

Following Fabrizio jij, Kishine [p^ and Bourbonnais pi], we investigate the aymptotic behaviours in/around the 
gapped phase by RG flow of the couplings. The partition function Z — J De^ where the action consists of one-particle 
hopping process 5'(i) and two-particle scattering 5(2) in the ladder, S ~ 5(i) -I- <S'(2). The dispersion is linearized with 
respect to the bonding(B) and anti-bonding(A) band which lie above(inner) the bonding band(LB, LA, RA, RB), 
where L, i? denote the right-moving, left-moving fields. Each band ranges by the bandwidth cutoff Kq. The interladder 
repulsion generates scattering vertices of Fig. 3. The action for the one-particle process is given 

^(1) ^Y.Y. \Gl\{K\^LXAK)UAK) + G-^^{K\^Rt,AK)Rr,o{K) (Al) 

K,c!ri=A,B 

where K = (fcy, A:j^,ze), K\\ — (fc||,ie) with the fermion frequency e^ = (2n -f l)7r/3, and the Green's functions are 
Glti = [i^-q + vp{k\\ + kpri)]^"^ , Gbi) — [i^-q — i'_F(fc|[ " fc_F)))]~"'^i with the Fermi momenta of the anti-bonding band 
kpA ~ kp ^ t±/vp and of the bonding band kps = kp + t±/vp. The action for the two-particle processes is given 



V,K,a 

+ gt{ai<J2a3<J4)L+,,{Ki)R+,AK2)L^as(.K3)R^,,AK4)S{Ki + K2 - K^ - K4) 

+ gt{<Ji<J2<yi<y4)L+,, {Ki)R+^, {K2)R^,,, {K3)L^,,, {K^)5{K^ +K2-K3- K^)] (A2) 

where the dimensionless coupHngs 5(fTicr20'3cr4) = g^^' S{ai(73)S{a2<J4) — g^^'S{ai(73)6{a2<J4) and go is for the intraband 
process, gf for interband forward scattering, gj, for interband backward scattering, gt for interband tunnehng. We scale 
the two-spin processes to remove the high-energy components, by parametrizing the bandwidth cutoff K(l) = K^e^^ 

with length I. The scaling equations for the scattering vertices are from — ^^di ~ ^ai^) ~ 22(1) where i — b for 
backward scattering and i — f for forward scattering. The vertex correction diagrams to the second order of the 
couplings(Z^(/)) are shown in Fig. 4. The linearized bands consist of the slow modes(s) (from the central portions of 
the branches) and the fast modes(f) (from the outer edges). The action is decomposed 'S'(i)_s + ^(2),s + ^(i)j + ^C2).f- 
Expanding the fast components exp[S'(i)j -I- >S'(2),/] with respect to <S'(2)./ and integrating out with the fast modes 
over dK{l) — ^K{l)dl, the partition function becomes Z — J_Dsexp[S'(i) .5 + S'(2),s + X];^o J\ ^^ ^{2) f >>c] where 
<<>>c denotes the connected diagram. The rcnormalized propagator and hopping amplitudes give renormalized 
action. 

G-,H^||) = [1 + mdl + 0(dP)]G-^i(K) (A3) 

gl^[l + zl,{l)dl + Oid\')]^^ (A4) 

il=[l + 0(dl2)]]ti (A5) 

where a stand for (o, f, b, t) and i for backward(l) and forward(2) scatterings. The Feynman diagrams contributing 
to the order dl are evaluated. Those in Fig. 4 give ^^(O? while the self-energy diagrams give z{l) of the scaling 
equations (Fig. 5). 

z(0 = .g« (0^ + g^^^ {If + gf' {ir + gf {if + gi'^ (if + 9^ H? + 9^' (if 

+ 9f\lf 9%il)9^^\l) - 9f\l)9f\l) - 9\f\l)9f\l) - 9'^t'\l)9?\l) (A6) 

The scaling equations follow, provided small Kq < Akp ; 

^5^(0 - -2i9i>'\ir + 9't'\l)9l'\l)) + 2{gf\l)gf\lf - gf\l)g^,'\l)gl'\l)) 

- 2g^\l){g^\lf + gf\lf + gf\lf + gf\lf - g^\l)gf\l)) (A7) 

^5^(0 = -(.9^(0^ +5r^(0^ + 5r^(0^) + 2(4^^ (03f^(0^ -5P^ (03^^(03^(0) 

+ 29^(1)9? H) - 9'^ - 9fil)9'o'\l) - 9r(l)9f{l) 

- 2g';^\l){gf\lf + gf\lf ~ g^\l)gf\l)) (A8) 

j/f\l) - ~2igf\ir + gi'\ir - gi'\l)gi'\l)) + 2[gl,'\l) g^Hlf g'^\l)9i'\l)9['\l)] 

^ 29f\lU'\ir+gf\ir+gi'\lfgi'\lf-9i'\l)gi'\l)] (A9) 

f/Pil) = -{9? (If 9?\lf) + 2{g^^\l)9?'(lf 9i'\l)9?\l)9l^\l)) 
+ 2g^l\lfgf\l) - gi'\lfg^,'\l) - gi'\lfgf\l) - ^(^^(O^ 
-2igf\l)igi'\lf+gf\lf-g['\l)gi'\l)) (AlO) 



y^\l) = -2{9'^\r)9f{l)+9fil)i'il)+i'il)9'^\l)-9'^\l)9'^\l) 
+ 9'iHl)9'i\l) - 3^(05^(0) - 2gi^\0[(5r(0 ~ 9f{l)f 

- 9'^\l)9^\l)+9'^\l)9f{l)+9?{l)9':^\l) - 9T^l)9fil)+9'^\l?+9Ti^m 

+ 2g'i\l){gf\lf-g'^\l)gf\l)+gf\lf) (All) 

j/:\l) - -2{9i'\l)9f\l) ^ 9^^\l)9^^\l) + 9^^\l)9f\l)) 

- 29i'\l)9i,'\l)gf\l) - 2gf )(/)[(.g^^)(/) - gf\l)r - g^'^(/).gf (/) 

+ .9^(04' '(0 + 9f\l)9o'\l) 9f\l)9f{l) + (5^'^(0 - 9f{l)f] 

- 2gi'\l)ig['\ir - gi'\l)g['\l) + g['\ir) (A12) 

^gf'(0 = -2(2g(i)(05(l)/(0-.9pH03f^(0-4'^(05f'(0 +5^(03^^(0 
+ .gP^(/).9f (0) + 2(2.g(i)(0ff^')(03f ^(0 - 9^'\l)9i'\l)9f\l) 
- 9f\l)9i'\l)9'^\l)) ^9^\l){9'^\l?+9^\l?+9f\lf+9f{l? 
+ 9^{lf+9f\lf-9^\l)9'i\l)~9f{l)9f\l)~9^{l)9f\l)) (A13) 

j,9'^\l) = -2{9i'\l)9['\l) - .9^ (0.9f '(0 - 3^(05^^(0) 

+ 2 (25^^) (04') (0.9^^ (0 - 4'^ (05^ (Offf ^ (0 - ff^'^ (Osf (Offf^ (0 ) 

+ 2gi'\l)gf\l)g['\l) - gi'\l)gf\l)gl^\l)) ~ 25f \/)(.9«(0^ + fffC/)^ + ^^^(0^ 

From the initial condition ga (0) = t^ /Jirvp = / > 0, the scahng equations give the strong-couphng fixed point ; 
^W ^ _;^^ ^(1) ^ Q^ ^(1) ^ -^^ ^(2) ^ 1/4(2/- 3), .9^^^ = 1/4(2/4-1), .9^ ^ = 1. However, the coupUngs actually do not 
flow to this fixed point, as far as the interchain hopping is greater than the gap size, since interchain pair fluctuations 
are dominant excitations Q . The scale invariant coupling for the charge part is gp^± = {g^ — 2g\^ ')±{g\' — 2gf ) . 
The invariant coupling for the spin part gcr,± = ffo ±5f where the + sign is for the total mode and — for the relative 
mode. From above g's, the spin stiffness is determined K^j^ = J y^. So, Ka,+ = means that the system developed 

the spin gap phase, while the charge part remains gapless from Kp^± — (Luther- Emery phase) ||2^. In our equation 
(8), the sine-Gordon mass term identifles W ^ [Jvp/t^)'^ , f^ , where g\ respresents some strong coupling value 

resulted from the scaling equations. 
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FIG. 1. spin ladder example : a two-leg ladder system. The dashed circles are for spin singlets (triplets, if excited) and for 
holon pairs after doped. 
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FIG. 2. phase diagram : the energy scale vs. the ratio W = [S±t± — Jx]/[5||t|| — Jy] 
"AF SL" stands for antiferromagnetic spin hquid which occurs ai W — Wo = J±/J\\ with 5 — 0. The spin gapsize A is 
marked for the point where gapped spinon phases branch from AF SL state, (a)gapped magnons; excitonic bound states 
of spinous, Luther-Emery liquid(LE) with gapless charge excitation (b) LE with gapful spinous and gapless holons, CDW (c) 
unstable solitons, gauge bosons, SDW (d) Luttinger liquid (paramagnetic); gapless excitations of both-charge and spin (e) 
superconducting (with d-wave symmetry) (f) electron phase: spinon and holon confined 




FIG. 3. interaction vertices. The solid lines are for the right-going fermions and the dashed lines for the left-going ones, g^^' 
and g^^' denote backward and forward scattering couplings, respectively. 
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FIG. 4. vertex diagrams to the next-to-the-leading order. 
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